Scaling theory for diffusion limited cluster aggregation in a porous medium 
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A scaling theory is developed for diffusion-limited cluster aggregation in a porous medium, where 
the primary particles and clusters stick irreversibly to the walls of the pore space as well as to each 
other. Three scaling regimes are predicted, connected by smooth crossovers. The first regime is at 
low primary particle concentrations where the primary particles stick individually to the walls. The 
second regime is at intermediate concentrations where clusters grow to a certain size, smaller than 
the pore size, then stick individually to the walls. The third regime is at high concentrations where 
the final state is a pore-space-filling network. 
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Diffusion-limited aggregation (DLA) and diffusion- 
limited cluster aggregation (DLCA) attracted much at- 
tention from the early 1980s following the discovery of 
scale- invariant (fractal) structures [l|, 0, 9- In DLA, 
primary particles undergo random walks before adher- 
ing irreversibly to a growing cluster Q. In DLCA, both 
primary particles and growing clusters undergo random 
walks, adhering to each other when they come into con- 
tact Despite a couple of decades or so of intense 
study, there still remain interesting questions about these 
processes. On the one hand deep puzzles linger about 
the origins of the scale invariance and whether a simple 
fractal picture suffices [4|, and on the other hand gener- 
alisations to more complicated situations have been con- 
sidered, such as DLCA as a function of primary particle 
concentration [f| or in the presence of sedimentation Q 
or flow @,||]. 

In this note I present a scaling analysis of DLCA as 
a function of the primary particle concentration for the 
case where the primary particles and clusters can adhere 
irreversibly to the walls of the porous medium as well as 
to each other. I consider the case where the primary par- 
ticles are much smaller than the pores. My motivation 
is to gain a deeper understanding of the deposition of 
small particles into porous materials, with applications 
that range from deep-bed filtration 9|, to the aggrega- 
tion and deposition of colloidal asphaltene in reservoir 
rocks in oil recovery (Toj |. The results presented here 
are also hopefully intrinsically interesting, although the 
model still lacks certain features present in real systems 
such as flow. This will be discussed further at the end. 
Previously, DLA has been studied in confined geometries 
fill ], but this is not really relevant to the present problem. 

Let me start by summarising what is known about 
DLCA [i2{ . Firstly, as already mentioned, the clusters 
have a scale invariant structure such that the aggregation 
number Af of a cluster of linear size R obeys 



N ~ (R/a ) d f 



(1) 



where ag is the primary particle size and df is the fractal 
dimension. This scaling law has been confirmed both 
in experiments on aggregating colloids and in computer 
simulations [H, [llj]. For the remainder of this note my 



focus will be on the practicallyrelevant case of three 
dimensions, where df as 1.75 (3l.[l3j. 

A second line of enquiry into DLCA has been into 
the kinetics of the process 0, [H, HI]. A long time 
ago, Smoluchowski presented a mean-field theory to de- 
scribe the evolution of the cluster size distribution [l7l ]. 
When updated to take account of the fractal structure 
of the aggregates, Smoluchowski's theory has turned out 
to be suprisingly accurate. One result of Smoluchowski's 
theory is that the cluster size distribution remains bell- 
shaped [Til Il5| . In scaling terms this justifies the no- 
tion of a characteristic cluster size. Armed with this 
foreknowledge, it is possible to use a simple scaling ar- 
gument to determine how the characteristic cluster size 
grows with time [l(|. In fact this is the key to the de- 
velopment, since the scaling argument generalises in a 
straightforward way to the problem of DLCA in a porous 
medium, by taking account of the depletion kinetics of 
clusters as they adhere to the walls. 

The kinetic scaling argument for pure DLCA runs as 
follows. Let the n(t) be the cluster number density at 
time t. Then, in mean field theory, dn/dt = — fc agg n 
where fc agg is the aggregation rate for the characteris- 
tic clusters. Again taking a mean field approach, this 
aggregation rate is given by the classic Smoluchowski 
diffusion-limited reaction rate fc agg ~ RDn, where R is 
the characteristic cluster size and D is the correspond- 
ing diffusion coefficient. This result can be derived by 
considering the steady state diffusive flux to a sphere of 
radius ~ R from a concentration field of number density 
ri at infinity. Now, the Stokes-Einstein relation implies 
D ~ (do / R) x Dq where D is the primary particle diffu- 
sion coefficient. One concludes that dn/dt ~ —aoDon 2 . 
Integrating this gives 1/n — 1/no ~ aoDgt, where hq is 
the number density of primary particles. Defining the 
aggregation number for the characteristic clusters to be 
N = nrj/ro, one finds 



N = 1 + k Q t 



(2) 



where kg ~ ao-Do n o is the characteristic aggregation time 
for primary particles. 

Remarkably, this is exactly the same answer as can be 
obtained from the full Smoluchowski theory. For k$t 3> 1 
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FIG. 1: Double-logarithmic plot showing scaling regimes for 
diffusion-limited cluster aggregation in a porous medium, here 
represented schematically as a slit. In order of increasing 
volume fraction, the regimes are primary particle deposition, 
cluster deposition, and formation of a pore-space-filling net- 
work. Parameters are the primary particle size ao, the char- 
acteristic pore size A > flo, the primary particle diffusion 
coefficient Do, and the fractal dimension df. In three dimen- 
sions df ~ 1.75 and the slopes of the lines are l/(d/ — 1) ~ 1.33 
and -3/(3 - df) « -2.40. 

it implies TV ~ t z with z ~ 1 being the kinetic scaling 
exponent. This has been confirmed both in experiment 
and in computer simulations [HI, [lj| • 

The fractal nature of the growing clusters means that 
they eventually expand to fill space, forming a network 
or gelled state. The time scale for this can be worked 
out by considering the effective cluster volume fraction 
4> ~ nR 3 . Assuming k t 3> 1 and using Eqs. §T§ and 
one finds <fi ~ ^(kot) 3 ^^ 1 where 4>q ~ na\ is the volume 
fraction of primary particles. Since df < 3, <p increases 
with time and gelation occurs when <j> — > 1. The gelation 
time is therefore 

V-^/DoX^'. (3) 

A point to note in idealised DLCA is that the final state of 
the system is always a gelled state. Experimentally other 
factors may intervene of course, such as sedimentation or 
creaming, or cluster rupture due to flow. 

The framework has now been set for the generalisation 
to DLCA in a porous medium. In this more complex 
situation, the basic idea is that there is a competition 
between the rate at which primary particles and clusters 
find each other, and the rate at which they are depleted 
by becoming stuck to the walls. I consider the case where 
the porous medium has a unimodal pore size distribution 
and therefore a characteristic pore size A 3> ao (the gen- 
eralisation to multiple length scales or fractal pore spaces 
is not too difficult). I also take the starting point to be a 
random dispersion of primary particles at a number den- 
sity no (volume fraction 4>q ~ rioOo). It might be objected 
that this is experimentally unfeasible but this point will 
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FIG. 2: Intermediate regime of Fig. [T] double-logarithmic 
plot showing how the cluster-cluster aggregation rate falls 
faster than the cluster-wall depletion rate as a function of 
time. Where the two rates cross over, clusters become at- 
tached to the walls, halting the aggregation process. The time 
at this cross over point should be less than the gelation time, 
otherwise wall cluster deposition is pre-empted by gelation. 
In three dimension, df ~ 1.75 and so 1/dj « 0.57. 

be discussed in more depth at the end. 

Since X 2 /D is the time it takes for the cluster of size 
R to diffuse a distance of order the pore size, the rate 
at which such clusters are depleted by sticking to the 
walls is given in scaling terms by fcd op ~ D/X 2 . This 
assumes R <C A, which will be verified a posteriori. This 
rate should be compared to the aggregation rate from 
Smoluchowski's theory above. Both are changing in time 
and to compare on a like-for-like basis the t-dependence 
should be explicitly extracted. One finds 

fcagg ~irb' ^c P ~^(i+M)- i/d '. (4) 

Two situations should be considered, depending on the 
initial values of the rates for primary particles. Firstly, if 
&agg ^ fcdop at t = 0, then the primary particles diffuse 
to the walls much faster than they undergo aggregation. 
This condition implies ko <C -Do/A 2 or 

0o « (ao/A) 2 . (5) 

In this limit aggregation is irrelevant and the final state is 
one where primary particles are stuck to the walls. More- 
over, since there are no particles per unit volume spread 
over the walls at a surface area per unit volume ~ 1/A, 
the area fraction covered by particles is ag?ioA ~ ^oA/ao- 
But, using Eq. ([5]), </>oA/a <C ao/A <C 1. In other words 
the primary particles are very sparsely distributed over 
the walls. Finally the time to reach the final state is 
Tppd ~ A 2 /Z?o where Dq is the primary particle diffusion 
coefficient. This can be written as 

Tppd ~ a 2 ,/ A) x (A/a ) 2 - (6) 

This is independent of the primary particle concentration 
of course. 
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Now consider the case where fc agg » fcd ep at t = 0. 
In this situation, the primary particles start to aggre- 
gate before they see the walls. Eq. (j4j shows that, as 
time progresses, fc agg falls faster than fc<jep, since df > 1. 
Eventually the wall depletion rate must overtake the ag- 
gregation rate. This is strongly suggestive of a crossover, 
in other words one expects aggregation to proceed as in 
standard DLCA until fc agg ~ &dcp, a ^ which point the 
clusters are removed by sticking to the walls. This is il- 
lustrated in Fig. [5] Making this ansatz in Eq. (H|) , and 
rearranging, shows that this point is reached on a time 
scale 

r cld ~ al/D x (A/a ) 2 VW-D . ( 7 ) 

Once the clusters are stuck to the walls, no further 
changes occur, so this is also the time to reach the final 
state of the system. The argument assumes fcoTcid ^ 1 
but one can easily show that this is true. In fact, one 
can show that r c id — > T pp d, and the aggregation number 
of the characteristic clusters N — > 1, as (f>o — > (cio/A) 2 
from above, implying there is a smooth crossover between 
the regime of primary particle deposition and the present 
cluster deposition regime. 

The time to reach the final state in Eq. ([7]) is an in- 
creasing function of 4>o and at some point it must surpass 
the gelation time in Eq. ([3|) which is a decreasing function 
of 4>q. If this happens, the system will form a pore-space- 
filling network before it reaches the state where clusters 
become attached to the walls. To avoid this fate requires 
that T c id <C Tgd , which after some rearrangement implies 
0o <SC (a /A) 3 ~ d/ . Thus the cluster deposition regime 
occurs only in an intermediate range of volume fractions, 

(a /A) 2 «0 O « (a /A) 3 -^. (8) 

Note that this intermediate regime exists since oq <C A, 
and df > 1 is required for connected clusters. 
Conversely, for 

(ao/A) 3 -* « 0o, (9) 

the final state is a pore-space-filling network or gel. The 
time to reach the final state in this situation is just the 
gelation time given by Eq. ([3]). 

At the crossover between cluster deposition and gela- 
tion, one has r c id ~ T gc i ~ o^/Dq x (A/ao) 3 . This is found 
by setting <f) ~ (a /A) 3 -^ in either of Eqs. © or (0). 

It is worth commenting on the properties of the clus- 
ters in cluster deposition regime of Eq. ([5]). After some 
algebra one finds that the terminal characteristic cluster 
size, in units of the pore size, obeys R/X ~ 1 / 2 where 
(j> ~ nR? is the effective cluster volume fraction. This 
shows that the clusters only grow to span the pores as 
one approaches the crossover to the gelation threshold 
where <$>—>!. This, incidentally, is also the promised 
a posteriori justification for R -C A mentioned above. 
It also eliminates a possible fate of the aggregating sys- 
tem: at no point does the characteristic cluster size grow 
large enough to span the pores (R ~ A), without gelling 



(cj) <C 1). Finally, by analogy to the primary particle de- 
position case, the effective area fraction covered by clus- 
ters in the cluster deposition regime is R 2 n\. But one 
has i? 2 nA - R 3 n(X/R) ~ - 1 / 2 . This means 

that clusters deposit individually in the cluster deposi- 
tion regime, and only start to touch each other as one 
approaches the gelation threshold. 

With these considerations it has now been demon- 
strated that the crossovers between all the regimes are 
continuous. At this point the analysis is complete. There 
are three predicted scaling regimes: a primary particle 
deposition regime at the lowest primary particle con- 
centration, a cluster deposition regime at intermediate 
concentrations, and a pore-spacc-filling network or gela- 
tion regime at the highest concentrations. This is shown 
schematically in Fig. [TJ using a slit to represent the pore 
space. In terms of the primary particle volume fraction, 
these regimes are given by Eqs. (0, ((HJ) and (J9j) respec- 
tively. Likewise the time taken to reach the final state 
in each of the regimes is given by Eqs. J6]), 10 and ([3]) 
respectively. A regime of pore-spanning clusters does not 
occur. 

Fig. [T] shows perhaps a surprising prediction, namely 
that the time to reach the final state increases with con- 
centration in the intermediate cluster deposition regime, 
reaching a maximum at the gelation threshold. Naively 
one might have expected T pp d to be an upper bound 
for the time to reach the final state. After all, non- 
aggregating particles would stick to the walls on a time 
scale ~ Tppd independent of the concentration. However, 
in the cluster deposition regime for an aggregating sys- 
tem, r P pd is initially much longer than the aggregation 
time. What happens then is that most of the primary 
particles form aggregates rather than finding the walls. 
The diffusion coefficient of an aggregate is smaller than 
that of a primary particle, hence the time to reach the 
final state (i. e. the time to diffuse to the walls) must 
exceed r pp d. 

The present note is a self-contained presentation of 
a scaling theory for DLCA in a porous medium. Al- 
though the predictions seem reasonable, and internally 
consistent, the approach is based at least in part on the 
ansatz, illustrated in Fig. [2l that there is a crossover 
between cluster-cluster aggregation and cluster-wall de- 
pletion. For future work it is obviously prudent to test 
this. An obvious way is to use computer simulations. 
For these, it is probably sufficient to consider DLCA in a 
slit of width A, as indicated schematically in Fig. [T] Ex- 
perimental tests could also be attempted although they 
are harder to make since a proper comparison should 
ideally start from randomly dispersed primary particles 
which are colloidally unstable. Usually, colloidal insta- 
bility is induced by adding salt, or by some other mix- 
ing process, and in a porous material it seems inevitable 
that this would involve flow fields. However, ingenious 
methods have been devised for the in situ destabilisation 
of colloidal suspensions, for example using the enzyme- 
catalysed hydrolysis of urea to increase the ionic strength 



4 



[TH . Another possible way might be to shine UV light on 
colloids which have been stabilised by photo-destructible 
surfactants [TH, although so far as I know this latter 
method has not been experimentally tested. 

Two other directions for future work are as follows. 
Firstly it would be interesting to include the effects of 
a reduced sticking probability, both for particle-particle 
collisions, and for particle-wall collisions. One might an- 
ticipate a variety of cross-overs, based on our understand- 
ing of how reaction-limited cluster aggregation (RLCA) 
crosses over to DLCA [l5j]. Secondly for realistic appli- 
cations it may be important to include the effects of flow 
fields. One can anticipate that flow would influence the 
behaviour in two distinct ways. Firstly, advection by 



the flow may produce a different scaling behaviour for 
clusters formed at Peclet numbers that are larger than 
unity (f|. Secondly, if the shear stresses are sufficiently 
large, they would have the effect of rupturing clusters 
and therefore limiting growth [8|. 
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